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Abstract. We studied the role of fundamental constants in an updated recombination sce- 
nario, focusing on the time variation of the fine structure constant a and the electron mass 
m e in the early Universe. Using CMB data including WMAP 5-yr release, and the 2dFGRS 
power spectrum, we put bounds on variations of these constants, when both constants are 
allowed to vary, and in the case that only one of them is variable. In particular, we have 
found that -0.019 < Aa/a < 0.017 (95% c.L), in our joint estimation of a and cosmo- 
logical parameters. Finally, we analyze how the constraints depend on the recombination 
scenario. 



1. Introduction 

Time variation of fundamental constants is a 
prediction of theories that attempt to unify the 
four interactions in nature. Many observational 
and experimental efforts have been made to 
put constraints on such variations. Cosmic mi- 
crowave background radiation (CMB) is one of 
the most powerful tools to study the early uni- 
verse and in particular, to put bounds on possi- 
ble variations of the fundamental constants be- 
tween early times and the present. 

Previous analysis of CMB data (earlier 
than the WMAP five-year release) includ- 
ing a possi ble variation o f a h a ve been per- 
formed bv iMartins et all d2002|) : iRocha et all 
(2003); llchikawa et all (I2006I): IStefanescJ 
( bo07|) : iMosquera et all d2008l) : [Landau et alj 
(2008) and including a possible variation of m e 
have been performed by llchikawa et al.1 (120061) : 



Scoccola et all ( 2008b): Landau et alJ (2008); 
Yoo&Scherrerld2003l) . 

In the last years, the process of recom- 
bination has been revised in great detail 
(Chlu ba & Sunvaevll2009albl: iHirata & Forbesl 
120091) . and in particular, helium recombination 
has been calculated very precisely, revealing 
the importance of considering new physical 
processes in th e calculation of the re c ombi- 
nation history dDubrovich & Grachevl 120051: 
Switzer & Hiratal l2008at IHirata & Switzerl 
2008t ISwitzer & Hiratal l2008bt 

Kholupenko et al.ll2007l) . ~ 

In a previous paper (IScoccola et al.l2 008a), 
we have analized the variation of a and 
m e in an improved recombination scenario. 
Moreover, we have put bounds on the possi- 
ble variation of these constants using CMB 
data and the power spectrum of the 2dFGRS. 
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In Section [2] we review the dependences on 
a and m e in the recombination scenario. In 
section [3] we present bounds on the possible 
variation of a a nd m P using the 5-yr data re- 
lease of WMAP dHinshaw et al.l2009h together 
with other CMB experi ments and t he power 
spectrum of the 2dFGRS (tCole et al.ll2005h . In 
additi on to the results of our previous work 
dScoccola et al.l l2008al) . here we also present 
bounds considering only variation of one fun- 
damental constant (a or m e ) alone. A compar- 
ison with similar analyses by other authors is 
presented in Section|4] 

2. The detailed recombination 
scenario 

The equations to solve the detailed recom- 
binati on scenario can b e found for exam- 
ple in IWong et ail d2008l) , as they are coded 
in Recfast. In the equation for helium re- 
combination, a term which accounts for the 
semi-forbidden transition 2 3 p-l's is added. 
Furthermore, the continuum opacity of HI is 
taken into account by a modification in the es- 
cape probability of the photons that excite he- 
lium atoms. 

The cosmological redshifting of a transi- 
tion line photon K and the Sobolev escape 
probability ps are related through the follow- 
ing equation (taking Hei as an example): 



1 



Hei 



gHeI,2'p « H eI,l's^2ip-li s -P S 



(1) 



where A HeI 2i p -ii s is the A Einstein coeffi- 
cient of the He I 2'p-l 's transition. To include 
the effect of the continuum opacity due to HI, 
ba sed on the approx i mate formula suggested 
bv lKholupenko et aU d2007l) . ps is replaced by 
the new escape probability p esc — p s + p m n,n 
with 
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Pcon.H 



1 + OHeT He 



where cr H ,i.,(vHei,2'p) is the H ionization cross- 
section at frequency VH e i,2'p an d Av Dj 2' P = 
Vnei^p V^^B^'M/'MHeC 2 is the thermal width of 
the He^p-l's line. 

The transition probability rates A HeI 2i p _ii s 
and A H eT73p_|i 5 can be expressed as follows 
(Drake & Morton||2007|): 



3c 2 



(4) 



where <Wy is the frequency of the transition, 
and i(J) refers to the init ial (final) state of the 
atom. It can be shown dScoccola et alj feOQSa) 
that to first order in perturbation theory, the de- 
pendence of the bracket goes as the Bohr ra- 
dius «o- On the other hand, <y, ; is proportional 
to the difference of energy levels and thus its 
dependence on the fundamental constants is 
coy - m e a 2 . Consequently, the dependence of 
the transition probabilities of Hei on a and m e 
can be expressed as 



A™ - m e a\ 



(5) 



The dependences on a and m e of all the phys- 
ical quantities relevant at recombination are 
summarized in Table [TJ 

Fig. [TJ shows the ionization history for dif- 
ferent values of a. Recombination occurs at 
higher redshift if a is larger. On the other hand, 
there is little change when considering differ- 
ent recombination scenarios, for a given value 
of a. Something similar happens when varying 
m e . 

With regards to the fitting parameters an e 
and hue, it is not possible to determine yet 
the effect that a variation of a or m e would 
have on them . Howe ver, we have shown in 
IScoccola et al.1 d2008al) that for the precision of 
WMAP data, there is no need to know these 
dependences. 



(2) 3. Results 



and 



We performed our statistical analysis by ex- 
ploring the parameter space with Monte Carlo 
Mark ov chains generated wi th the CosmoMC 

code dLewis & Bridld |2002|) which us es the 

87r 3/2 cr H>ls (y H ei,2ip)v 2 eI21 Av Dj2 i p (l - x p ) Boltzmann code CAMB dLewis et al.1 12000) 



W P 2'p-l' 



(/He - *HeIl)c 2 
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Description 


Physical Quantity 


Dependence 


Binding Energy of Hydrogen 




a 2 m e 


Transition frequencies 


nt2s>V H eI,2 1 8' V HeI^ 3 s 


(t 2 m f 


Photoionization cross section n 


(T n (Z,kv) 


-1 -2 

a m e 


Thomson scattering cross section 


(Tj 


2 -2 

a m e 


Recombination Coefficients Case B 




3 -3/2 

am,. 


Ionization Coefficients 


y0H,y6 H eI 


a 3 


Cosmological redshift of photons 




a~ 6 m~ 3 


Einstein A Coefficients 


'-' 


a 5 m e 


Decay rate 2s — > Is 


A H , A H ei 


a & m e 



Table 1. Dependence on a and m e of the physical quantities relevant during recombination. 



new recombination scenario 
standard recombination scenario 




Fig. 1. Ionization history allowing a to vary with time. From left to right, the values of ^- are 
1.05, 1.00, and 0.95, respectively. The dotted lines correspond to the standard recombination 
scenario, and the solid lines correspond to the updated one. 



and Recfast to compute the CMB power spec- 
tra. We modified them in order to include the 
possible variation of a and m e at recombina- 
tion. Results are shown in Table |2] 

We use data from the WMAP 5 -year tem- 
perature a nd temperature-p olarization power 
spectrum dNolta et al.1 120091) . and other CMB 
exper iments such as CBI (Readhead et al. 
120041) . ACBAR ( |Kuoetal.l |2004 . and 
BOOMERANG Jpiacentini et all l2006t 



I Jones et al . 2006), and the powe r spectrum 
of the 2dFGRS dCole et al l 120051) . We have 
considered a spatially-flat cosmological model 
with adiabatic density fluctuations, and the 
following parameters: 



P = |n B /i 2 ,Q C DM/? 2 ,©,T,— , 



Aor Am„ 



■,n s ,AA (6) 



where Qcdm^ 2 is the dark matter density in 
units of the critical density, © gives the ratio 
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of the comoving sound horizon at decoupling 
to the angular diameter distance to the surface 
of last scattering, t is the reionization optical 
depth, n , the scalar spectral index and A s is the 
amplitude of the density fluctuations. 

In Table [2] we show the results of our 
statistical analysis, and compa re them with 
the on es we have presented in lLandau et al.l 
(2008), which were obtained in the stan- 
dard reco mbination scenario ( i.e. the one de- 
scribed in (SeageretaDllOOO), which we de- 
note PS}j_jffldjisingWMAP3 dffinshaw et al.1 
120071: iPage et al.l 120071) data. The constraints 
are tighter in the current analysis, which is 
an expectable fact since we are working with 
more accurate data from WMAP. The bounds 
obtained are consistent with null variation, for 
both a and m e , but in the present analysis, 
the 68% confidence limits on the variation of 
both constants have changed. In the case of a, 
the present limit is more consistent with null 
variation than the previous one, while in the 
case of m e the single parameters limits have 
moved toward lower values. To study the ori- 
gin of this difference, we perform another sta- 
tistical analysis, namely the analysis of the 
standard recombination scenario (PS) together 
with WMAP5 data and the other CMB data 
sets and the 2dFGRS power spectrum. The re- 
sults are also shown in Table [2] We see that the 
change in the obtained results is due to the new 
WMAP data set, and not to the new recombi- 
nation scenario. In Fig.|2]we compare the prob- 
ability distribution for Aa/ao in different sce- 
narios and with different data sets. In Fig. [3] 
we do the same for Am e /(m e )o- Bounds on the 
fundamental constants are shifted to a region of 
the parameter space closer to that of null varia- 
tion in the case of a. On the other hand, limits 
on the variation of m e are shifted to negative 
values, but still consistent with null variation. 
The bound on Q/,/z 2 is also shifted to higher 
values. 

We present here the results of our statistical 
analysis when only one fundamental constant 
is allowed to vary together with a set of cos- 
mological parameters. We obtained these re- 
sults using data from WMAP5, CBI, ACBAR, 
BOOMERANG, and the P(k) of the 2dFGRS. 
The constraints (with 1-cr errors) on the vari- 



ation of a are Aa/ao = -0.002 + 0.009 in the 
standard recombination scenario, and Aa/ao = 
-0.001 + 0.009 in the detailed recombination 
scenario. For m e , both bounds are Am e /(m e )o — 
-0.01 + 0.03. The limits are more stringent 
than in the case of joint variation of the con- 
stants. This is to be expected since the parame- 
ter space has higher dimension in the later case. 
The values for the cosmological parameters are 
consistent with those from the joint variation 
analysis. 



4. Discussion 

The obtained results for the cosmological 
parameters are in agreement within lcr with 
the ones obtain e d by th e WMAP collaboration 
dDunklev et al.l 120081) . without considering 
variation of fundamental constants. It is 
also interesting to compare our results with 
the works of iNaka shima et al. (2008]) and 
Menegoni et al.1 ( 20091) where only the varia- 
tion of a was analized using WMAP 5 year 
release and the Hubble Space Telescope (HST) 
prior on the Hq- In these works, the HST prior 
on Ho is used to reduce the large degeneration 
between Hq and a and find more stringent 
constraints on a variation . Howe ver, we have 
shown in Mosquera et al.l (120081) (where the 
variation of a was analysed using the WMAP 
3 year release), that more stringent constraints 
on a can be found using the power spectrum 
of the 2dFGRS. Indeed, our constraints on a 
al one are more stringent than those reported 
by lNakashima et al.l d2008l) . On the other hand, 
our constraints are of the same order tha n 
those presented by Menegoni et al ] (l2009h . 
using data from higher multipoles reported 
by recent CMB experiments, such as QUAD 



(OUaD collaboration: M. L. Brown et a 



2009) and BICEP dChiang et al.l 120091) . 
Finally, it is important to stress, that when 
using the HST prior on Ho, the correct value 
to be used is the value obtained using only the 
closest objects, since bounds obtained using 
other objects could be affected by a possible a 
variation. 
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parameter 


wmap5 + NS 


wmap5 + PS 


wmap3 + PS 


Q b h 2 


"■"'"-0.1X1)14 


"•"^^ -0.00085 


0218+ , .!!!!! 




1070 +0 ,S 

" -± " '"-0.0078 


"•*"' -0.0080 


o 106+ , ." 

"• ""-0.011 





" J -0.023 


1 03261+ .- ,? 4 

" J " -0.023 


l 033 + ?™ 

-0.029 


T 


0870 +0 ,™" 

"•""'"-0.0081 


" "°" -0.0084 


" " v "-0.014 




"•"" -0.015 


003 + !!°.^ 

"•"" J -0.015 


-0 023+ .- , 2 ? 

"•"'-'-0.025 


Affl„/(m,)n 

e i v c so 


-0.019 +0 ,°^ 

-0.049 


-0.017 + o^; 

-0.051 


0.036+ , ,™ 

-0.078 


n s 


962 +00 ' 4 


963 +0 015 


" y '"-0.019 


A s 


q r\CQ+0.042 
J -O.IMI 


3 05203+ 004269 


3 054 +0 073 

Jl -1.073 


H 




70 3 +(U 

'"•-'-6.0 


70.4^ 



Table 2. Mean values and l<x errors for the parameters including a and m e variations. NS stands 
for the new recombination scenario, and PS stands for the previous one. 




A a / Oq A a / a 

Fig. 2. One dimensional likelihood for — . Left figure: for WMAP5 data and two different recom- 

*> ff b 

bination scenarios. Right figure: comparison for the standard recombination scenario, between 
the WMAP3 and WMAP5 data sets. 
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